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ABSTRACT 



Non-local interactions are assumed for the deuteron with form factors 

9c(k) = jo(b\k) jofak) for the central part responsible for the S-state and 

gr(k) = ji(bik) jifok) for the tensor part responsible for the D-state, where b\ 
and 62 are range parameters for the proton and neutron, respectively. The analytically 
obtained wavefunctions in coordinate space have different forms in three different 
regions. The inner most region is between r = and r = 62 — 61 (assuming 
b2 > 61), followed by a region between r = 62 — b\ and 61 + 62 and finally the region 
r > 61+62- The resulting wavefunctions and their derivatives are found to be continuous 
at the boundaries. The ensuing calculations are simplified by setting b\ = 62 = b. 
Good agreement is obtained for the quadrupole moment. Neutron-proton scattering 
calculations will follow in the next communication. 
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1. Introduction 



Non-local interactions for the neutron-proton scattering and bound state (the 
deuteron) are written in the form 

V c (k,k') = -^L gc (k)gc(k% (1.1) 

and 

V T (k,k') = ^-g T (k)g T (k>), (1.2) 

where Vc(k, k') is the central part responsible for the S-State of the deuteron and 
Vr(k,k') is the tensor part corresponding to the D-state of the deuteron. Ac and Ay 
are the strengths of the central and tensor parts of the potential, respectively and are 
assumed positive. The form factors gc(k) and gr(k) are assumed to be of the form 



9c{k)=j Q {bik)j {b 2 k), (1.3) 

and 

g T (k)=ji(b 1 k)j 1 (b 2 k), (1.4) 

where ji(bk) is a spherical Bessel function of order /, b\ and 62 are range parameters for 
the proton and neutron, respectively. The deuteron (reduced) S- and D- wavefunctions 
in momentum space, u(k) and w(k), respectively are then [1,2] 



4) = i/MW)») 

V n k 2 + a 2 



and 
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where A and B are the normalisation constants for the S- and D- wavefunctions, 
respectively and a 2 = M\E^\, M is mass of the nucleon, is the energy of the 
deuteron. 

Using 

oo 

u(r) = y^r J k 2 dku(k)j (kr), (1.7) 
o 

oo 

w(r) = y — r J k 2 dkw(k)j2(kr) : (1-8) 
o 

and assuming 62 > 61, 

(i) for r < b2 — b\ 

u(r) = Aio(abi) ko(ab2) (ar)io(ar), (1-9) 

u;(r) = -Bii(o;6i) ^1(0:62) (or ^(ar), (1-10) 

(ii) for 62 — &i < y < &i + &2 

u (r) = 9 A , {1 - cosh[a(& 2 - bi)]e~ ar - e -^ bl+b ^ sinh(ar)}, (1.11) 
2a z b\b2 

B r ~ 2/7,2 , .2\ , 3 4/,2 i2\2 , ,.2/72 , t2\ oi , 2 2 



W(r) = 16o^2 { 2 « W + b 2) ~ 4 - ^2 [« ( & 2 - + M&i + $ - 8] + a 2 r 
+ 8ar k2(ar)[(a 2 bib2 — 1) cosha(&2 — 61) + a(^2 — &i) sinha(&2 — &i)] 
- 8ari 2 («r)e" Q(&1+b2) [a 2 6i6 2 + a(&i + 62) + 1]} 



(1.12) 



(iii) for r > 61 + 62 

u(r) = Aio(abi) 20(0:62) (or) ko(ar), 

w(r) = B i\(abi) i\(ab2) (ar) k2(ar). 
The modified spherical Bessel functions ii(x) and fc/(x) are given by 



sinh a; 

io{x) = , 

x 



and 



x cosh x — sinh x 

ii(x) = 2 ; 

x z 



. , . (x 2 + 3) sinhx — 3x cosh a; 

12 (X) = ^ , 



h(x) = 



X 



e x (x + l) 
k i( x ) = =2 ■ 



. , , e~ x (x 2 + 3x + 3) 

k < X ) = Z3 ' 



To simplify the calculations using these wavefunctions, set b\ 
coordinate space wavefunctions then become 



(i) for < r < 2b 

A 



u(r) = — ^ {1 - e~ ar - e- 2ab sinh (ar)} (1.21) 
2a z b z 



(1.22) 



W(r) = 16^F { - 4(1 - " V) + " V + (1 - ^ 

- 8(1 - a 2 b 2 ){ar) k 2 (ar) - 8(1 + abf e~ 2Q& (ar) i 2 (ar)} 

(ii) for r > 2b 

u(r) = Ail(ab)(ar)k (ar), (1.23) 

w(r) = Bi\{ab) (ar)k 2 (ar), (1.24) 
The S-state and D-state probabilities are calculated using 

oo 

Ps = J k 2 dk\u(k)\ 2 , (1.25) 
o 

oo 

P D = J k 2 dk\w{k)\ 2 , (1.26) 
o 

with the result 



A 2 

P S = [8a& - 9 + 4e- 2ab (2ab + 3) - e^ ab {Aab + 3)], (1.27) 

= [56« 5 6 5 - 135a 4 6 4 - 80a 3 6 3 + 450« 2 6 2 - 315 

- 60(ab + l) 2 (2a 3 6 3 + 3a V - 7) e - 2ab - 15(ab + l) 3 (4aV + 7ab + 7)e" 4a6 ]. 

(1.28) 

The normalisation constants A and B are, of course, related by the relation 



P s + P D = 1. (1-29) 

Initial estimates, based on fitting the quadrupole moment and the root mean 
square radius show that B 2 3.0 A 2 for a range parameter of b = 1.475 fm. Using 
a = 0.23165 fm^ 1 , then ab = 0.342. This results in S-wave and D-wave probabilities 
of approximately 96% and 4%, respectively. The resulting values for the normalisation 
constants are A = 0.905 fm -1 / 2 and B = 1.57 fm -1 / 2 to 3 significant figures. The 
asymptotic normalisations are then 

A s = Ail(ab), (1.30) 

and 

A D = Bi\(ab), (1.31) 

resulting in the values A s = 0.941 fm _1 /2 an( j a d = 0.0208 fm _1 /2. The D/S ratio is 
then rj = Ad /As = 0.022. Figures 1 and 2 show u(r) and w(r) as compared with the 
NIJM I potential model calculations [3]. 

The root mean square radius is defined by 



_ 1 

frms — ^ 

and the quadrupole moment by 



r 2 dr[\u(r)\ 2 + \w(r)\ 2 ], (1.32) 



Q = — J r 2 w(r) [V8u(r) — w(r)]dr. (1.33) 
o 

Using Gaussian quadratures, one can easily perform the integrals numerically to obtain 
r rms = 2.08 fm and Q = 0.286 fm 2 . 



2. Conclusions 



Spherical Bessel functions were used in non-local interactions describing the 
deuteron. It was shown that the resulting wavefunctions exhibit different behaviour in 
three different regions of coordinate space. This property is due to using the spherical 
Bessel functions. The agreement with the quadrupole moment of the deuteron was to 
2 significant figures. More accuracy is expected when b\ is set to be different from hi- 
Scattering calculations and the 3 Si phase shifts are deferred to the next communication. 
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